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A method for calculating the theoretical heat consumption in glassmaking has been developed for the zonal 
method of calculating exterior heat transfer. Computational results are presented for heat consumption in a 
280 - 420 tons/day container-glass furnace. 
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In solving the exterior problem of heat transfer in 
gas-electric glass furnaces the source term (- Q gf ) in the 
nonlinear equations of heat transfer and heat balance for the 
surface zones of the molten glass with boundary conditions 
of the second kind is calculated from the expression [1] 

8gf = Ql.l + 8i.2 + 0m.mt~ Qe’ (1) 

where Q x , and Q x 2 are, respectively, the heat consumption 
in glass formation and heating of the molten glass, kW; Q m mt 
represents the heat losses through the melting tank masonry, 
kW; and, Q e is the additional electric heating (AEH), kW. 

To calculate heat transfer by the resolvent zonal method 
[2] the equation determining Q x x [5] is put into the form 

8\,\^ Pf c G b (ci^ + q Xi ), ( 2 ) 

where P fc is the furnace capacity, kg/sec; G h is the batch con¬ 
sumption per 1 kg molten glass, kg/kg; q z is the total heat of 
the glass-forming reactions per 1 kg batch, kJ/kg. 

Here 

9ij = ?i.i + 9u + 9i.3 + 9i.4» 

where q x x , q x 2 , q x 3 and q x 4 are, respectively, the heat con¬ 
sumption in evaporation of moisture and heating the pri¬ 
mary melt and products of degasification and glass melting, 
kJ/kg. 

1 First President of Russia B. N. El’tsin Ural Federal University, 
Ekaterinburg, Russia. 

2 E-mail: vdzuzer@yandex.ru. 


The batch consumption per 1 kg molten glass is 


where G is the yield of molten glass per 1 kg batch, 
kg/kg. 

In turn 

G mg =10-0.01G p . d +G' C 
and 


c 100-G" ’ 

where G p d is the total yield of the products of degasification 
per 1 kg batch, kg; G' c is the cullet consumption per 1 kg 
batch, kg/kg; G" is the cullet fraction in the total mass of 
batch and cullet, %. 

The total heat of the glass formation reactions per 1 kg 
batch is determined from the relation 

<7i =0.0ljG,4„ 

i=i 

where G ( is the mass of the 7th oxide, kg; q t is the specific 
heat consumption in the reaction producing the 7th oxide, 
kJ/kg; and, n is the number of oxides. 
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The mass of the zth oxide per 100 kg batch is found from 
the expression 


maximum temperature of the molten glass in the melting 
tank 


G,-=i)G b . m * 0 (l0-H 2 O), 

;=i 

where G b is the content of the z'th material in the raw materi¬ 
als used in the batch, %; K 0 is the content of the main sub¬ 
stance or corresponding oxide of the z'th component, frac¬ 
tional content; ffO represents the moisture content of the 
batch, fractional content; and, k is the number of batch com¬ 
ponents containing the zth oxide. 

The heat consumption in the evaporation of moisture is 

q x .\ =^h 2 o, 

where K = 2514 is the heat consumption in the evaporation 
of 1 kg moisture, kJ/kg. 

The heat consumption in heating primary melt is 
q 1.2 ^'mg ^mg f 

where c mg is the specific heat capacity of the primary melt at 
Z‘ mg , kJ/(kg • K); t mg is the maximum heating temperature of 
the molten glass in the melting tank, °C; and, t b is the batch 
temperature at the time of loading into the furnace, °C. We 
shall take t mg = 1500°C and t b = 40°C. 

The temperature dependence of the specific heat of the 
primary melt is given by the equation 

c mg =4.187x (0.1605+ Llx 10 4 ? mg ). 

The heat consumption in heating the products of degasi¬ 
fication is given by the relation 

z/i.3 = 0-01kp d c pd t pd , 

where V d , c pd , and t pd are, respectively, the volume (per 
100 kg batch), m 3 , the average specific heat, kJ/(m 3 • K), and 
degasification temperature, °C. We shall take ? pd = 1500°C. 

The volume of the products of degasification is calcu¬ 
lated from the equation 


Qi. 2 =P ic G c c c (t mg -t c ), (3) 

where G c is the cullet consumption per 1 kg molten glass, 
kg/kg; c c is the specific heat of the cullet, kJ/(kg • K); and, 
Z c is the temperature of the cullet at the time of loading into 
the furnace, °C. We take t c = 40°C. 

The specific heat of the cullet is calculated from the re¬ 
lation 

c c =4.187x (0.1794 +0.632x 10 4 ? mg ). 

The cullet consumption per 1 kg molten glass is 


or 

G c = G b G' c . (4) 

Substituting the expression (4) into Eq. (3) we obtain 

01.2 = 'Pfc^b ( Jc c c( ? mg _? c )• (5) 

Summing the left- and right-hand sides of Eqs. (2), (5) 
and rearranging we obtain an expression for calculating the 
theoretical heat consumption in glassmaking Q t , kW, 

0 theor =' P p.d ( Jb[(9x + <7l.l + ?1.2 +#1.3 + <7l.4) + 

Gc c c0mg ~t c )l ( 6 ) 

We shall now put Eq. (6) into a form suitable for the 
zonal method of calculating the external heat transfer: 

0theor = P p.d(4mg + ?h)> ^ 

where q mg is the specific consumption of heat on glass for¬ 
mation, kJ/kg, and q h is the specific heat consumption in 
heating the primary melt and cullet, kJ/kg. 

In turn, 


where G ; is the total yield of the zth component of degasifica¬ 
tion kg; p ( is the density of the zth component of degasifica¬ 
tion, kg/m 3 . 

The heat consumption in glass melting is 
q l 4 = 347.52x (L0-0.01G p d ). 

In the heat balance for the melting part of the furnace the 
quantity Q x 2 represents the heat carried away by the working 
flow of the molten glass [4]. In the zonal calculation of heat 
transfer Q ] 2 is the heat consumption in heating cullet to the 


0mg +< 7l.l +9l,3 +( ?1.4)i (8) 

9h =G b [q l2 +G' c c c (t mg -f c )]. (9) 

Separating the consumption in glass formation (8) and 
heating (9) in the total theoretical heat consumption in glass¬ 
making (7) makes it possible differentiate their distribution 
along the melting tank [3]. 

To illustrate the method of determining g theor we shall 
examine the computational results obtained for the heat con¬ 
sumption in glassmaking for a 280 - 420 tons/day con¬ 
tainer-glass furnace (melting tank area 115.77 m 2 ). The 
chemical composition of the glass is as follows (wt.%): 
71.9Si0 2 , 2.0Al 2 O 3 , 0.2Fe 2 O 3 ; 9.6CaO; 2.2MgO; 13.1Na 2 0; 
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TABLE 1. Batch Recipe for 100 kg Molten Glass 


Content by weight 

K.aw material 

Parts 

% 

Quartz sand, GOST 22551-71, PB-150 brand 

65.80 

55.995 

Nepheline - feldspar material, 



TU 5726-036-001193861, PShS-0.2-21 brand 

8.55 

7.276 

Ground dolomite, GOST 23672-79, 



DM-20-0.1 brand 

10.77 

9.165 

Chalk, technical, TU 113-08-667-98, A brand, 



grade 1 

10.37 

8.825 

Soda ash, technical, GOST 5100-85, A brand. 



highest grade 

20.09 

17.096 

Sodium sulfate, technical, GOST 6318-77, 



highest grade 

1.57 

1.336 

Chromium colorant, 



TU 2362-001-70470322-2006, A brand 

0.25 

0.213 

Charcoal, granular, TU 12.36.210-91 

0.11 

0.094 

Total: 

117.51 

100 


0.55K 2 O; 0.25Cr 2 O 3 ; and, 0.2SO 3 . The batch recipe is pre¬ 
sented in Table 1. 

The content of recycled cullet in the total batch and 
cullet mass is 30%. The moisture content of the batch is 4%. 

Computational Results for the Intermediate Data 
for Eq. (6) as well as <jr mg and q b 


(7 p d- kg/100 kg batch.19.79 

Fp.d, kg/100 kg batch.12.68 

G' c , kg/kg batch.0.429 

G mg , kg/kg batch.1.231 

Gb, kg/kg glass mass.0.812 

G c , kg/kg molten glass.0.348 

c mg , kj/(kg • K). 1.3229 

c c , kj/(kg ■ K).1.1481 

Cp.d, kJ/(m 3 ■ K).2.1518 

</: ,. kJ./kg.100.56 

</.;■ k.l/kg. 2377.60 

<7i. 3 , kj/kg. 409.27 

</i 4 ■ k.l/kg. 278.74 

</>.. k.l/kg. 355.00 

9m g , kJ/kg. 928.58 

</,. k.l/kg.2514.52 


The computational results for the theoretical heat con¬ 
sumption in glassmaking (Table 2) in combination with the 
coefficients of heat transfer through the furnace masonry [5] 
make it possible to complete the initial database for deter¬ 
mining the source term (- Q mg ). 

We now revisit the data in Table 2. As the glass melting 
output increases, the theoretical heat consumption in 


TABLE 2. Computational Results for the Theoretical Heat Con¬ 
sumption in Glassmaking and Furnace Thermal Power 


Furnace 

capacity, 

tons/day 

Specific glass 
melting output, 
tons/(m 2 ■ day) 

£?theor ’ 

kW 

f?p.d > 
kW 

£?p.d £?theor ’ 

kW 

280 

2.42 

11158.2 

15006.9 

3848.7 

300 

2.59 

11955.2 

15590.3 

3635.1 

320 

2.76 

12752.2 

16049.9 

3297.7 

340 

2.94 

13549.2 

16607.4 

3058.1 

360 

3.11 

14346.2 

17156.2 

2810.0 

380 

3.28 

15143.3 

17703.8 

2560.5 

400 

3.46 

15940.3 

18257.3 

2317.0 

420 

3.63 

16737.3 

18824.0 

2086.8 


glassmaking increases as does the thermal power of the fur¬ 
nace, given by Q pd = BQ^ W +Q e , kW (B is the gas con¬ 
sumption in the furnace, m 3 /sec; glow is the lower working 
heat-generating capacity of the gas, kJ/m 3 ). In view of the 
known limitation on the temperature of the dinas roof of the 
furnace (1580 - 1600 °C) the entire increment to the thermal 
power cannot be realized because of flame heating. It is dis¬ 
tributed between the heat of combustion of the fuel (BQ ^ W ) 
and the power of the electric-heating of the melt in the melt¬ 
ing tank ( g e ). Determining the efficient value of AEH for 
different values of the specific glass melting output is a topi¬ 
cal problem of mathematical modeling of glass furnaces. 
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